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Abstract 

We investigate a mechanism for cooling a lead based on a process that replaces hot electrons by cold ones. The central 
idea is that a double quantum dot ■with an inhomogeneous Zeeman splitting acts as energy hlter for the transported 
electrons. The setup is such that hot electrons with spin up are removed, while cold electrons with spin down are added. 
The required non-equilibrium condition is provided by the capacitive coupling of one quantum dot to the shot noise of 
a strongly biased quantum point contact in the tunnelling limit. Special attention is paid to the identification of an 
operating regime in which the net electrical current vanishes. 
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1. Introduction 

The program of devising novel refrigeration schemes for 
electronic systems defines a variety of fundamental prob¬ 
lems of potential technological interest. The pumping of 
heat from a cold to a hot electron reservoir by suitably ma¬ 
nipulating the electron transport may be viewed as a form 
of rectihed motion [niiE]. Electron pumps can be adia¬ 
batic msiii] or (if they use inelastic transitions) nonadia- 
batic [71 [H]. The latter case represents a particular instance 
of ac-driven electron heat transport mm- It has been 
recognized that the effect of nonadiabatic driving can be 
very similar to that of an energetic, far-from-equilibrium 
energy source that promotes inelastic transmission of elec¬ 
trons through a spatially asymmetric setup [imiiiiiiii]. 
In this paper, we follow the approach of these works and 
investigate the possible use of the shot noise of a nearby 
quantum point contact (QPC) as the nonequilibrium en¬ 
ergy source that, through capacitive coupling, induces rec¬ 
tification in a spatially asymmetric device which in our 
case will be a double quantum dot (DQD) system. 

In mesoscopic physics, recent experimental progress in 
energy harvesting m turned the fundamental issue of 
heat balance into a topic of practical interest [16) . In 
particular, quantum dot setups have been argued to pro¬ 
vide a convenient framework for the controlled transport 
of heat between electron reservoirs [3 m HD HH], as has 
been experimentally confirmed |19j . On the other hand, 
it has also been argued that capacitive coupling to sys¬ 
tems with fluctuating charges can be a source of energy 
and of induced heat transport [5D], a claim that has been 
recently observed |21j . A similar prediction for quantum 
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cavity systems |22] has also received experimental confir¬ 
mation [231 El] • Analogous physics has been predicted for 
small electron systems interacting through the exchange 
of phonons [25) or microwave photons [26) 

When considering electron exchange between reservoirs, 
one typically assumes each lead in a grand canonical state 
determined by its chemical potential and its temperature. 
Then thermal excitations involve occupied states above 
the Fermi energy and unoccupied states below. Our focus 
lies on the associated excitation energies. Consequently, 
we speak of cooling when electrons in states above the 
Fermi energy are removed or when holes below are filled. 
Thus, one way to cool a lead is to contact it with another 
lead at a different chemical potential while applying an 
appropriate energy filter which can be realized, e.g., by 
the gap of a superconductor m or by a resonant level. 
The same principle can be applied to the cooling of of 
a two-dimensional electron gas mediating the charge flow 
between two leads, as has been predicted |28| and observed 
[29) ; for a review see [30) . A recent experiment is presented 
in Ref. [31) . which is based on the mechanism proposed in 
Ref. EOj- 

The above cooling scheme, however, relies on an applied 
bias which also causes a net electron transfer. In particular 
when one lead is actually a small grain, one soon reaches a 
situation in which the grain becomes electrically charged 
and the cooling process comes to an end. In this article 
we follow Ref. [7] by proposing a mesoscopic heat pump 
that avoids a net charge transport while operating between 
two leads with equal chemical potentials. In consistency 
with the second law of thermodynamics, such a heat flow 
from cold to hot requires some non-equilibrium condition. 
In Ref. |7], this has been theoretically achieved by an ac 
gating. Here by contrast, we employ the shot noise of a 
nearby quantum point contact. The shot noise of a nearby 
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conductor has already been investigated as a source of rec¬ 
tification [311 [331133] . Here we adapt the study in Ref. 
|3i] . where the fluctuating conductors is a QPC in the 
tunnel limit, to the pumping of heat with zero net electric 
current between the cold and the hot leads. The advan¬ 
tage of using a charge-fluctuating QPC is two-fold: First, 
its broad-band excitation is less sensitive to small detun¬ 
ings. Second, most recent quantum dots already include 
such a QPC (33 ES] so that our proposal can be realized 
readily. 

The behavior of quantum dot systems subject to the 
effect of nearby fluctuations can be studied as a particu¬ 
lar case of the “particle-bath” problem, where attention 
is paid to the effect on a few physical variables (or on a 
reduced Hilbert space) of many degrees of freedom that 
are mathematically traced out m- The effect of the dis¬ 
sipative environment is that of rendering the dynamics ef¬ 
fectively irreversible despite the time-reversal symmetry of 
the underlying microscopic dynamics. On a macroscopic 
level, irreversibility becomes manifest in the fact that heat 
can flow spontaneously only from a hot to a cold reservoir. 
Such a quantum dissipation approach underlies the mas¬ 
ter equation calculation described in section 2, where the 
leads degrees of freedom are traced out. 

This paper is arranged as follows. In section 2 we de¬ 
scribe the model considered and the calculation method 
employed. Section 3 contains analytical and numerical re¬ 
sults for the case of noninteracting electrons. In section 
4, we discuss the effect of Coulomb repulsion at the dots, 
focusing on the displacement of the particle-hole symme¬ 
try point. Section 5 is devoted to conclusions. Appendix 
A provides a detailed account of the arguments leading 
to the analytical results obtained for the charge and heat 
currents which are presented in section 3. 

2. Model and master equation 

2.1. DQD coupled to a QPC in the tunnel regime 

Our model sketched in Fig. consists of two capaci- 
tively coupled electric circuits, namely an unbiased DQD 
and a strongly biased tunnel contact. The latter entails 
non-equilibrium noise on the former and, thus, drives the 
DQD out of equilibrium. The DQD is modeled by the 
Hamiltonian 

Hdqu = ^ XI (cLc2<t -f cl^cia) + Y X 

CT—(l) 

2 — 1,2 K 

where i = 1,2 and tr = f, 4- label the two quantum dots and 
the spin degree of freedom, respectively, while the multi 
index k = {i, a) combines both. The terms describe spin- 
independent tunneling H, Coulomb repulsion U, an inho¬ 
mogeneous Zeeman splitting > 0, and a global shift 
of the onsite energies caused by a background gate volt¬ 
age (with the sign convention that a positive Igate shifts 




Figure 1: Sketch of the double quantum dot (DQD) coupled to a 
QPC in the tunnel regime. We focus on the heat balance in the right 
lead of the DQD, where hot spin-up electrons are removed while cold 
spin-down electrons are added. The coupling to the shot noise of the 
QPC creates a non-equilibrium situation that breaks the symmetry 
between the forward and the backward process. In the absence of 
the QPC, detailed balance inhibits cooling. 


the DQD levels upwards). For simplicity, we restrict our¬ 
selves to undetuned dot levels (in the absence of the mag¬ 
netic field) and assume that both intra-dot and inter-dot 
Coulomb repulsion are equal. 

Dot 1 is tunnel coupled to the left lead L described by 
.^dot-lead ^i.cr^h.q.cr) and dot 2 ac¬ 

cordingly to the right lead. The effective dot-lead coupling 
is given by the rate F = 2Tr J2q | 14 ,, 9 p^(e —Cg) which we as¬ 
sume independent of the energy e and the spin projection 
a. Then the equilibrium distribution of the lead electrons 
is given by the Fermi function fi{e) = [exp{e/kBTi) + 1]“^ 
where d = L,R, i.e., while setting p, = 0 for both leads, we 
allow for different lead temperatures 

Our second system is a QPC in the tunnel limit be¬ 
tween two leads modeled by the Hamiltonian iJqpc = 
Sfc ^kc\ck + J2k' ^k'c\,Cki, where k and k' label the modes 
of the left and the right lead, respectively, including the 
spin. The leads are weakly coupled by the tunnel Hamil¬ 
tonian A = A_|_ -b A_, where 

A+ — ^ ^ CyCk, (2) 

k,k' 

transfers an electron from the left to the right QPC 
lead, while A_ = A)(_ describes the opposite process. In 
the continuum limit, the matrix elements t^k' are en¬ 
compassed by the energy-independent QPC conductance 
G = 2,TrJ2kk' \tkk'\'^S{e — ek)S{e — Ck') in units of the con¬ 
ductance quantum Gg = e^/h. The electrons on dot 1 
enhance, via Coulomb repulsion, the barrier between the 
QPC leads and thereby reduce the tunnel matrix elements 
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tk^i. This effect is captured by a prefactor x = {1 — sni) 
in the tunnel Hamiltonian such that 

^QPC = + ^-) (3) 

accounts for both the QPC and its coupling to the DQD. 
The latter is mediated by the occupation of dot 1, ni = 
X^o-^icr- For consistency, the dimensionless coupling s 
must obey 0 < s < 1/2. 

2.2. Master equation 

Our theoretical description is based on the formal elimi¬ 
nation of all four leads such that we remain with a reduced 
master equation for the DQD. While the treatment of the 
leads coupled to the DQD follows a standard procedure, 
the systematic elimination of the tunnel contact is less 
common and has been performed only recently |38L 134) . 
In those works, a counting variable for the tunnel contact 
allowed to compute the full counting statistics and corre¬ 
lation functions, while here it is sufficient to compute the 
action of the QPC on the DQD. We sketch here the deriva¬ 
tion of the formalism of Ref. [351IM] as required for our 
present purposes. 

We start from the Liouville-von Neumann equation for 
the full density operator which we transform to the inter¬ 
action picture with respect to i?DQD and the lead Hamil¬ 
tonians. The remaining terms are treated within second- 
order perturbation theory in the dot-lead tunnelings and 
in the QPC tunneling to obtain the Bloch-Redfield master 
equation [39] 


Symmetrizing the time integral amounts to neglecting the 
energy renormalization stemming from principal values. A 
main ingredient is the correlation function of the QPC tun¬ 
nel operator, C{t) = (A(t)A(O)) = C+{t) + C-{t), where 
C±{t) = (Azp(t)A±(0)) is readily evaluated from its defini¬ 
tion and the assumption that the leads are voltage biased. 
In Fourier representation it reads [42] 

C±(w) i_exp[-(;ia;±eoF)/fcBT]' 

We restrict ourselves to the limit in which the QPC bias 
V is much larger than any other relevant frequency scale 
of the DQD. Then C{oj) = he^GV becomes frequency in¬ 
dependent so that C(t) oc 5{t). Then we obtain for the 
QPC Liouvillian the Lindblad form 

C.CiPCP = l{xpx-[x^,p\l2) (7) 

with the effective rate 7 = 27r/Qpc/eo. For a treatment 
beyond the large-bias limit, see Ref. [53] . 

In our numerical implementation of the master equa¬ 
tion formalism, we use the many-particle eigenstates of the 
DQD Hamiltonian as a basis and keep all off-diagonal el¬ 
ements of the density matrix. This ensures to capture the 
level repulsion stemming from the tunnel coupling, which 
is rather relevant for levels close to the Fermi energy. 

3. Heat balance in the absence of interaction 


P 




/ dt tneads[Vn,[Vn{-t),pf^ Rq]], 

n p 


(4) 


for the reduced DQD density operator p. Rq refers to the 
grand canonical ensemble of each lead, while the operators 
Vn represent the two tunnel contributions in 

f^DQD-ieads- Here, Vn(t) stands for Vn in the interaction 
picture. 

The evaluation of the Liouvillian /logo-leads for the in¬ 
coherent DQD-lead tunneling is rather standard, see e.g. 
the appendix of Ref. [ID] . It yields operators that describe 
jumps between many-particle DQD states differing by one 
electron. The transition rates contain Fermi functions re¬ 
flecting the initial occupation of the lead modes. To obtain 
an expression for the heat balance, we multiply in the su¬ 
peroperator for the electric current each tunnel process by 
the energy with respect to the chemical potential which 
the electron carries to the lead miiiiiiiin]. 

To obtain the Liouvillian for the action of the QPC, 
we evaluate the t-integral in Eq. 0 and find the QPC 
Liouvillian 


1 

>CQpcP=^y dt C{t)[x{-t)px + xpx{t) 


(5) 


—xx{—t)p — px{t)x\. 


To explain the central idea of the cooling mechanism, 
we consider the channels for the spin-up electrons and for 
the spin-down electrons separately. To do so, we ignore 
Coulomb repulsion which couples the spin-up and the spin- 
down channels. Later we will see that the interaction term 
in the Hamiltonian Q shifts the working point and affects 
the efficiency. 

3.1. Pumping mechanism 

We consider the DQD sketched in the upper half of 
Fig. [3 Let us focus on the channel for the spin-down 
electrons and assume equal temperatures, = Tj^. Sup¬ 
pose that an electron enters dot 1 from the left lead at 
the corresponding onsite energy. The electron may pro¬ 
ceed to dot 2 and eventually to the right lead. Due to 
the lower occupation at higher initial energy, the opposite 
process occurs with lower probability. Thus we expect a 
net transport that in the right lead fills holes below the 
Fermi surface, which corresponds to cooling. 

However, this picture is incomplete. Owing to the tunnel 
coupling between the dots, DQD eigenstates are delocal¬ 
ized, but in such a way that (spin-down) electrons in the 
ground state are mostly in the left dot while the excited 
ones dwell mostly in the right dot. Therefore, electrons in 
the excited states are more likely to leave to the right lead, 
while electrons preferentially enter the ground state from 
the left lead. The quantitative analysis [see Eqs. ^ and 


3 



(§ below] reveals that, when both leads have the same 
temperature and chemical potential, those processes fail 
to yield a net effect, i.e., both the electric current and the 
heat current vanish. This is of course what one expects for 
such an equilibrium situation. In more technical terms, the 
full description obeys detailed balance, as it should. 

To obtain any net transport, we must drive the system 
out of equilibrium. For this purpose, we couple the DQD 
to a strongly biased QPC in the tunnel regime; see the 
Hamiltonian ([^. Then the shot noise of the tunnel cur¬ 
rent induces transitions between the ground state and the 
excited state of the DQD. Since the ground state is more 
strongly populated, we witness a net excitation and, conse¬ 
quently, the transport process from left to right dominates 
the one in the opposite direction. We thus observe both 
charge pumping and cooling of the right lead. 

The channel for the spin-up electrons behaves similarly. 
It yields a net transport of hot electrons from the right 
to the left lead, which again corresponds to cooling. In 
the symmetric situation I4ate = 0 sketched in Fig. the 
transported charges of both spin channels compensate each 
other, while the energy transfer adds up. The net result is 
cooling of the right lead without charge accumulation. 


3.2. Quantitative analysis of the individual channels 
In order to substantiate the above discussion, we have 
solved the master equation for a single spin channel within 
the rotating-wave approximation. This approximation 
means that, within a Pauli-type master equation approach, 
off-diagonal density matrix elements are neglected and 
only the populations of the eigenstates are considered. 
The QPC-induced transitions between the (single-particle) 
ground state and the excited state of the DQD have been 
treated as a perturbation. Moreover, we have ignored the 
doubly occupied state. While details of the calculation 
can be found in [Appendix A( we present here the results 
for the electric current and the heat balance in the limit 
7 F: 


1 fe-h 

4 1-/g/e (e2+f72)3/2’ 


^ _ 7 /e - /g 

8 1-/g/e(e2+D2)l/2 


+ (b^ate 


( 8 ) 

(9) 


where e is the Zeeman energy gradient, e = {B 2 — Bi)f2, 
and the Fermi functions at the QDQ energies are abbrevi¬ 
ated as fg = f{Eg - fi) and fe = f{Ee - yd). 

If the ground state lies well below the Fermi surface, so 
that /g ~ 1, while /e is significantly smaller, the prefactor 
involving the Fermi functions becomes unity. Then both 
the electric current and the heat balance depend only on 
the DQD configuration and the distance to the symmetry 
point. 

An important implication of Eqs. § and 0 is that both 
vanish in the absence of the QPC, since then 7 = 0 . This 
underlines that our treatment is consistent with detailed 
balance. 



Figure 2: (a) Current and (b) heat balance as a function of the 

gate voltage Fgate for various DQD temperatures and zero Coulomb 
interaction. Parameters are F = 2£7 = B\/10 = B 2 I 6 = = 

eoV/80, QPC conductance G = 1, and coupling strength s = 0.1. 
The dotted lines mark the analytical solution given in |Appendix A| 
for kgT = F/lO. Inset: Zoom of the heat balance in the region close 
to Fgate = 0 where Q < 0. 


3.3. Cooling 

Our next goal is to provide numerical results for the dif¬ 
ferent operating regimes. In doing so, we plot in Fig. [^the 
electric current and the heat balance as a function of the 
gate voltage and for parameters that otherwise correspond 
to the sketch in Fig. [l] 

For very negative gate voltage, all four levels lie below 
the Fermi surface and no transport occurs. When shifting 
all levels upwards, the excited state of the spin-up channel 
will cross the Fermi function and a pump current from 
the right to the left lead sets in. The initial energy of the 
corresponding electrons in the right lead is far below the 
Fermi surface. Therefore the creation of holes corresponds 
to heating the right lead, which visible in the large positive 
value of Q at Igate = — 5F. For less negative Vgate, the 
energy of the created holes is also less negative and Q 
is diminished. Once the ground state also lies above the 
Fermi energy, the process comes to rest. For positive Vgate, 
the spin-down channel is active and adds hot electrons. 
While this leads to an opposite electric current, it also 
corresponds to heating. 

For zero temperature, the borders between the different 
regions are sharp and both the current and the heat bal¬ 
ance between the peaks vanish exactly. For larger temper¬ 
atures, the curves smear out. While at Igate = 0, the elec¬ 
tric current must vanish for symmetry reasons, the heat 
balance assumes negative values (see inset of Fig. [^. This 
means that we find a cooling process in which hot electrons 
of the right lead are replaced by cold ones. 
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Figure 3: Heat current as a function of the temperature of the left 
lead for (a) equal temperatures T = Tl = Tr and (b) for pumping 
against a temperature gradient for various Tr. The Zeeman splitting 
in dot 2 is B 2 = hF, while all other parameters are as in Fig.[^ The 
arrow in panel (a) marks the minimum of the curve for Bi = 12r; 
the vertical lines in panel (b) mark the zeros of the heat current and, 
therewith, the temperature up to which cooling can be achieved. 






Figure 4: (a) Heat production rate on the right lead as a func¬ 

tion of the Zeeman splittings B\ and B 2 at the symmetry point for 
T = 3r /ks and other parameters as in Fig. The cross line indi¬ 
cates its minimum for positive Zeeman splittings. The solid contour 
line indicates vanishing heat current, (b) Corresponding slices at 
constant B 2 = lOF (solid line), B 2 = 20F (dashed), and B 2 = 30r 
(dashed-dotted). 


Figure l^a) shows the heat balance as a function of the 
temperature for different Zeeman splittings. Interestingly, 
at zero temperature, Q — Q irrespective of the magnetic 
field gradient, i.e., not only the cooling but also the heating 
vanishes as expected. Moreover, the data indicate that 
cooling is possible only on the side on which the splitting 
is smaller. In the symmetric situation (dashed curve for 
Bi = B 2 ), we already observe significant heating in the 
right lead. For symmetry reasons it equals the heating 
of the left lead. Therefore, the coupling to the out-of¬ 
equilibrium QPC augments the total thermal energy, a 
behavior that agrees with our expectations. 

A more systematic investigation of the magnetic field 
dependence is provided in Fig. The overview as a func¬ 
tion of Bi and B 2 indicates that the optimal cooling of 
the right is found for B 2 clearly smaller than Bi, but not 
too small. One might have expected that Bi should be 
as large as possible to ensure a strong population of the 
relevant states in the left lead. However, in such a case 
the onsite levels are strongly detuned so that the effective 
transition matrix element of the perturbation between the 
DQD eigenstates becomes small; see [Appendix A[ 

A most intriguing question is whether one can reduce 
the thermal energy of one lead even when it is already at 
a lower temperature than the other lead, i.e., whether one 
can pump heat from cold to hot. In Fig. |^b) we show 
the results for pumping heat against a temperature gradi¬ 
ent. Globally, we find that this is possible for moderate 
temperature differences of roughly 10%. 


4. Coulomb repulsion 


To make our study applicable to realistic quantum dots, 
we have to include Coulomb repulsion. For 17 = 0, we 
had chosen our working point at Vgate = 0 where both 
channels are symmetric with respect to the Fermi energy. 
Therefore we start by writing the interaction term in the 
DQD Hamiltonian 0 in a more symmetric form with the 
help of the identity 


k^k' 


= E 



+ (V-l)y]n„-JiV(iV-l), 

K 


( 10 ) 


where Uk = cJIck is the occupation of the single particle 
level K = 1,... ,7V. The interpretation of this identity is 
that our interaction Hamiltonian can be expressed by the 
particle-hole symmetric terms cJIck — 1/2 = — (c^c^ ~ 1/2) 
plus an onsite energy. To be specific, the interaction term 
in Eq. 0 is particle-hole symmetric if a gate voltage shifts 
the TV = 4 levels by —317/2 [recall the prefactor U /2 in the 
interaction term of the Hamiltonian 0]. 

The predicted shift of the operating point is indeed visi¬ 
ble in the current and the heat balance as a function of the 
gate voltage and the interaction strength plotted in Fig.[^ 
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Figure 5: (a) Electric current and (b) heat balance as a function 

of the gate voltage Vgate and the Coulomb interaction U for k^T = 
r/2. Other parameters are as in Fig.j^ Black contour lines indicate 
vanishing current and heat current, respectively. The grey dashed 
line, Vgate = —3C//2, in panel (b) marks the points for which the 
electric current is zero. That line tracks the particle-hole symmetry 
point which is displaced by the interaction. 



Figure 6: (a) Current and (b) heat current as function of the gate 

voltage Vgate as in Fig. but for nonzero Coulomb interaction U = 
5r. Vertical lines indicate the transition energies listed in [Appendix] 
[A| Inset: Heat current in the region close to the symmetry point 
Vgate = —(3/2)f/ which is indicated by a vertical line. 



u [r] u [r] 


Figure 7: (a) Heat current on the right lead at the symmetry point, 

as a function of the DQD temperature and the Coulomb interaction 
for parameters as in Fig. The cross line indicates its minimum, (b) 
Corresponding slices at constant T = ST/kg (solid line), T = 2TI kg 
(dashed), and T = T/kg (dashed-dotted). 

The plot also reveals that with increasing interaction, each 
current peak splits into two peaks with a distance U. The 
values of these shifts can be appreciated in the horizontal 
slices of both panels shown in Fig. The physical rea¬ 
son for the peak separation is that for each channel, the 
onsite energies are relevant for the empty channel, while 
for the occupation with a further electron, one must over¬ 
come the Coulomb repulsion. This effectively augments 
the excitations energies by U. 

As a drawback, the interaction energy not only shifts 
the working point, but also changes Q quantitatively. The 
data in Fig. shows that the cooling effect has a maxi¬ 
mum at finite U where the location of the optimum value 
depends on the temperature. Thus, small interactions may 
be beneficial for cooling, but with increasing U, the cooling 
power decays and may even turn into heating. Thus an ex¬ 
perimental realization of our scheme should be attempted 
with quantum dots that have a large capacity. 

5. Conclusions 

We have shown that a double quantum dot subject to 
an asymmetric Zeeman splitting and capacitively coupled 
to the charge fluctuations of a nearby quantum point con¬ 
tact in the tunnel limit can experience rectification. In a 
suitable region of parameter space, the system can act as a 
heat pump that cools the lead coupled to the dot with the 
weaker Zeeman splitting, even if that lead is already the 
coldest one. Such a refrigeration scheme can operate with 
a zero net electric current across the double quantum dot 
system. This is possible thanks to a particle-hole symme¬ 
try that can be preserved if the gate voltage is suitably ad¬ 
justed to compensate the effect of Coulomb repulsion. For 
this setup to act as a heat pump, electric current creating 
shot noise must flow through the nearby quantum point 
contact. However, that electric current does not need to 
run parallel to the line connecting the two heat reservoirs 
through the double quantum dot system. This feature can 
become an important advantage in potential applications. 
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I = y'F cos(20)- 


a — b 


- ab + {2r + a + b)-f' 
s^TC{-S) b-a 

4 -ab+(2T+ a + 6 ) 7 ' 6^ 


, (A.7) 


Appendix A. Analytical solution for an individual 
channel 

Transport through an individual spin channel can be 
described by a Pauli-type rate equation P = A4P for the 
occupation probabilities P = {Pq, Pg, Pe)^ in the energy 
basis |0)(0|, | 5 )(( 7 |, |e)(e|. The ground and excited states are 
parametrized [23] by 

Ig) = — sin0|l)-|-cos6*|2), |e) = cos6*|l)-|-sin6*|2), (A.l) 


and 

<5= 

s 2 rC'(-^) b-a 

“ 8 -ab+(2T+ a + b)-f' ~T~ 

+ (^ate ~ (-^-8) 

The ground energy Eg and the excited energy occurring 
in the Fermi functions are given by Pg/e = Pgate — 
with 


where |1) and |2) refer to the left and right DQD state, 
respectively, in the local basis. The geometrical factors 
are determined by cos( 20 ) = and sin( 20 ) = |n|/i5, 

with level splitting <5 = |np and detuning 

^ sgn(cr), (A.2) 

where sgn(tT) is positive for cr = f and negative, otherwise. 
The Liouvillian reads 


M = 


—a — b 
a 
b 


a 

i — 

7' 


b 


j' 

-b-i 


(A.3) 


where a = F — a = r/(Pg — /r) and b = T — b — Tf{Ee — fi) 
are finite dot-lead tunneling rates. Albeit the leads are in 
equilibrium, electrons can still be excited by the driving 
of the QPC, which is modeled by Fermi’s golden rule rate 
7 ' = (s/ 2 )^ sin^( 20 )C'(—5) of the coupling operator Pqpc- 
From the current operator 


= -kU ± ( 5/2 _ -^1 + -^2 ( a . 9 ) 

The latter includes a displacement by a multiple k of the 
Coulomb interaction. To emphasize explicitly the depen¬ 
dence on the spin a and this multiple fc, we write in the 
following the current of the individual channels as / = Ik,a-- 
The two-channel case can be approximated by summa¬ 
tion of the single-channel solutions for spin-up and spin- 
down. For vanishing Coulomb interaction, the two channel 
current is obtained from J 07 -I- / 07 , while it is composed 
of /07 -I- /17 -I- / 2 ,t + for finite U S. The analytical 
heat current for the two channels case can be analogously 
defined. The two channels case for vanishing Coulomb 
interaction is in Fig. compared to the numerical solu¬ 
tion of the full master equation. Further, figure shows 
the transition energies for finite Coulomb interaction— 
the vertical lines, from left to right, are given by the en¬ 
ergies at (/c,cr,s) = (3,t,-),(3,t,+),(2,t,-),(2,t 
) +)> (1) —)) (1)4-) +); (0, —), ( 0 , 4 ,, -I-). 


J = 


0 asin^(0) bcos^{9) 

-asi-n^iO) 0 0 

—bcos^{9) 0 0 


and the heat current operator 


(A.4) 


Jq 


—asin^( 0 ) 
6 cos^( 0 ) 


asin^( 0 ) 

0 

0 


—5cos^(6*) 
0 
0 


(A.5) 


follows directly the current I = tr ^7P®* and the heat cur¬ 
rent Q — trj7QP®*, respectively. Hereby, tr = (1,1,1) 
denotes the trace operator and 


pst 


1 

-ab+i2T+ a + 6 ) 7 ' 


'ah -f (a -f b)Y 
ab+{a + b)i 
ab+{a + b)i 


(A. 6 ) 


References 

[1] M. Biittiker, Transport as a consequence of state-dependent dif¬ 
fusion, Z. Phys. B Condens. Matter 68 (1987) 161. 

[2] P. Reimann, Brownian motors: Noisy transport far from equi¬ 
librium, Phys. Rep. 361 (2002) 57. 

[3] P. Hanggi, F. Marchesoni, Artificial brownian motors: Control¬ 
ling transport on the nanoscale, Rev. Mod. Phys. 81 (2009) 387. 

[4] T. E. Humphrey, R. Newbury, R. P. Taylor, H. Linke, Reversible 
quantum brownian heat engines for electrons, Phys. Rev. Lett. 
89 (2002) 116801. 

[5] J. E. Avron, A. Elgart, G. M. Graf, L. Sadun, Transport and 
dissipation in quantum pumps, J. Stat. Phys. 116 (2004) 425. 

[6] J. P. Pekola, F. Giazotto, O.-P. Saira, Radio-frequency single¬ 
electron refrigerator, Phys. Rev. Lett. 98 (2007) 037201. 

[7] M. Rey, M. Strass, S. Kohler, F. Sols, P. Hanggi, Nonadiabatic 
electron heat pump, Phys. Rev. B 76 (2007) 085337. 

[8] F. Chi, Y. Dubi, Microwave-mediated heat transport in a quan¬ 
tum dot attached to leads, J. Phys.: Condens. Matter 24 (2012) 
145301. 

[9] M. Moskalets, M. Biittiker, Floquet scattering theory of quan¬ 
tum pumps, Phys. Rev. B 66 (2002) 205320. 


7 


















[10] L. Arrachea, M. Moskalets, L. Martin-Moreno, Heat produc¬ 
tion and energy balance in nanoscale engines driven by time- 
dependent fields, Phys. Rev. B 75 (2007) 245420. 

[11] J. P. Pekola, F. W. J. Hekking, Normal-metal-superconductor 
tunnel junction as a brownian refrigerator, Phys. Rev. Lett. 98 
(2007) 210604. 

[12] D. Segal, Single mode heat rectifier: Controlling energy flow 
between electronic conductors, Phys. Rev. Lett. 100 (2008) 
105901. 

[13] F. Sols, Aspects of quantum cooling in electron and atom sys¬ 
tems, Physica E 42 (2010) 466. 

[14] B. Cleuren, B. Rutten, C. Van den Broeck, Cooling by heating: 
Refrigeration powered by photons, Phys. Rev. Lett. 108 (2012) 
120603. 

[15] H. B. Radousky, H. Liang, Energy harvesting: an integrated 
view of materials, devices and applications. Nanotechnology 23 
(2012) 502001. 

[16] B. Sothmann, R. Sanchez, A. N. Jordan, Thermoelectric en¬ 
ergy harvesting with quantum dots. Nanotechnology 26 (2015) 
032001. 

[17] R. Sanchez, B. Sothmann, A. N. Jordan, M. Biittiker, Correla¬ 
tions of heat and charge currents in quantum-dot thermoelectric 
engines, New J. Phys. 15 (2013) 125001. 

[18] D. Venturelli, R. Fazio, V. Giovannetti, Minimal self-contained 
quantum refrigeration machine based on four quantum dots, 
Phys. Rev. Lett. 110 (2013) 256801. 

[19] H. Thierschmann, M. Henke, J. Knorr, L. Maier, C. Heyn, 
W. Hansen, H. Buhmann, L. W. Molenkamp, Diffusion ther¬ 
mopower of a serial double quantum dot. New J. Phys. 15 (2013) 
123010. 

[20] R. Sanchez, M. Biittiker, Optimal energy quanta to current con¬ 
version, Phys. Rev. B 83 (2011) 085428. 

[21] H. Thierschmann, R. Sanchez, B. Sothmann, F. Arnold, 

C. Heyn, W. Hansen, H. Buhmann, L. W. Molenkamp, Three- 
terminal energy harvester with coupled quantum dots. Nature 
Nanotech (2015) in press. 

[22] B. Sothmann, R. Sanchez, A. N. Jordan, M. Biittiker, Rectifi¬ 
cation of thermal fluctuations in a chaotic cavity heat engine, 
Phys. Rev. B 85 (2012) 205301. 

[23] B. Roche, P. Roulleau, T. Jullien, Y. Jompol, I. Farrer, 

D. Ritchie, D. Glattli, Harvesting dissipated energy with a 
mesoscopic ratchet. Nature Commun. 6 (2015) 6738. 

[24] F. Hartmann, P. Pfeffer, S. Hofling, M. Kamp, L. Worschech, 
Voltage fluctuation to current converter with coulomb-coupled 
quantum dots, Phys. Rev. Lett. 114 (2015) 146805. 

[25] O. Entin-Wohlman, Y. Imry, A. Aharony, Three-terminal ther¬ 
moelectric transport through a molecular junction, Phys. Rev. 
B 82 (2010) 115314. 

[26] T. Ruokola, T. Ojanen, Theory of single-electron heat en¬ 
gines coupled to electromagnetic environments, Phys. Rev. B 
86 (2012) 035454. 

[27] J. Pekola, Low-temperature physics: Tunnelling into the chill. 
Nature 435 (2005) 889. 

[28] H. L. Edwards, Q. Niu, G. A. Georgakis, A. L. de Lozanne, 
Cryogenic cooling using tunneling structures with sharp energy 
features, Phys. Rev. B 52 (1995) 5714. 

[29] J. R. Prance, C. G. Smith, J. P. Griffiths, S. J. Chorley, D. An¬ 
derson, G. A. C. Jones, I. Farrer, D. A. Ritchie, Electronic re¬ 
frigeration of a two-dimensional electron gas, Phys. Rev. Lett. 
102 (2009) 146602. 

[30] F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, J. P. 
Pekola, Opportunities for mesoscopics in thermometry and re¬ 
frigeration: Physics and applications. Rev. Mod. Phys. 78 
(2006) 217. 

[31] J. V. Koski, A. Kutvonen, T. Ala-Nissila, J. P. Pekola, On-chip 
Maxwell’s demon as an information-powered refrigerator (2015) 
arXiv:1507.00530 

[32] A. Levchenko, A. Kamenev, Coulomb drag in quantum circuits, 
Phys. Rev. Lett. 101 (2008) 216806. 

[33] R. Sanchez, R. Lopez, D. Sanchez, M. Biittiker, Mesoscopic 
coulomb drag, broken detailed balance, and fluctuation rela¬ 


tions, Phys. Rev. Lett. 104 (2010) 076801. 

[34] R. Hussein, S. Kohler, Capacitively coupled nano conductors: 
Ratchet currents and exchange fluctuation relations, Ann. Phys. 
(Berlin) (2015) in press; doi:10.1002/andp.201500141. 

[35] T. Ihn, S. Gustavsson, U. Gasser, B. Kiing, T. Muller, 
R. Schleser, M. Sigrist, 1. Shorubalko, R. Leturcq, K. Ensslin, 
Quantum dots investigated with charge detection techniques. 
Solid State Commun. 149 (2009) 1419. 

[36] D. Taubert, D. Schuh, W. Wegscheider, S. Ludwig, Determina¬ 
tion of energy scales in few-electron double quantum dots. Rev. 
Sci. Instrum. 82 (2011) 123905. 

[37] A. O. Caldeira, A. L. Leggett, Quantum tunnelling in a dissi¬ 
pative system, Ann. Phys. (N.Y.) 149 (1983) 374. 

[38] R. Hussein, J. Gomez-Garcfa, S. Kohler, Monitoring quantum 
transport: Backaction and measurement correlations, Phys. 
Rev. B 90 (2014) 155424. 

[39] A. G. Redfield, On the theory of relaxation processes, IBM J. 
Res. Develop. 1 (1957) 19. 

[40] R. Hussein, S. Kohler, Coherent quantum ratchets driven by 
tunnel oscillations: Fluctuations and correlations, Phys. Rev. B 
86 (2012) 115452. 

[41] U. Sivan, Y. Imry, Multichannel landauer formula for thermo¬ 
electric transport with application to thermopower near the mo¬ 
bility edge, Phys. Rev. B 33 (1986) 551. 

[42] G.-L. Ingold, Yu. V. Nazarov, Charge Tunneling Rates in Ul¬ 
trasmall Junctions, Vol. 294 of NATO ASI Series B, Plenum, 
New York, 1992, p. 21. 

[43] T. Brandes, Coherent and collective quantum optical effects in 
mesoscopic systems, Phys. Rep. 408 (2005) 315. 


8 


